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Abstract 

The process eq — > eq + 7 exhibits radiation zeros, i.e. configurations of the final- 
state particles for which the scattering amplitude vanishes. We study these zeros for 
both e + u and e + d scattering. The latter exhibits a type of zero which to our knowledge 
has not previously been identified. The observability of radiation zeros at HERA is 
discussed. 



1 Introduction 



In certain high-energy scattering processes involving the emission of one or more photons, 
the scattering amplitude vanishes for particular configurations of the final-state particles. 
Such configurations are known as radiation zeros. In the context of high-energy scattering, 
they were first discussed by Mikaelian, Sahdev and Samuel Q. To measure the magnetic 
moment fj,w of the W boson, they proposed and studied the W boson production processes 
ud — > W +, y and du — ► W'j. They found that the matrix element vanishes at a particular 
value of the c.m.s. frame scattering angle cos 9w = [&d(d) ~ e u(u)]/[ e d(d) + e u(u)] = — 1/3, 
independent of the photon energy A similar effect is seen in the W — ► qqj decay process . 
Experimentally, these radiation zeros have been observed recently by the CDF collaboration 
at the Tevatron pp collider. There has also been renewed theoretical interest, including 
studies on double photon emission processes pp — > W ±, yy — > i^wyy [ffl and the uniqueness 
of radiation zeros to the Standard Model @. The energy dependence of radiation zeros 
in pp — > 7 + X @ and using radiation zeros to probe the colour-charge of partons f| has 
also been studied. A review of recent developments in the subject can be found in Ref. |§. 

A first understanding of the phenomenon of radiation amplitude zeros was achieved in the 
pioneering work of Ref. ]nj . The vanishing of the scattering amplitude can be understood 
as arising from complete destructive interference of the classical radiation patterns of the 
incoming and outgoing charged lines in relativistic n-particle collisions. Taking the single 
emission of a photon as the paradigm process, it was shown that the amplitudes can vanish 
if the other particles participating in the process have the same sign of charge ej.[] 

Same-sign charge scattering occurs naturally in high-energy hadron collisions in subpro- 
cesses such as ud — > W +, y. However similar phenomena can be expected in lepton-hadron 
collisions, and in particular at HERA in processes such as eq — > eq + 7 for eq = e + u or e~d. 
Studies of radiation zeros for these processes at HERA were first performed by Bilchak [11 



and more recently by Doncheski and Halzen [JT2 



In a recent paper [|13[] we studied the distributions of soft gluon and photon radiation in 
eq — > eq scattering at HERA. The motivation was to demonstrate that the radiation patterns 
are different depending on whether the scattering takes place via standard t-channel 7*, Z* 
exchange or via the production of a new heavy, charged, colour-triplet 'leptoquark' (LQ) 
resonance in the s channel. Leptoquark production is one of several possible explanations 
for the apparent excess of high-Q 2 deep inelastic scattering events at HERA |14|, jl5| . 

A by-product of this study was the identification of radiation zeros in both the Standard 
Model and leptoquark e + q — ► e + q + 7 (q — u,d) scattering amplitudes. For a long-lived 
resonance (Flq — > 0) we found radiation zeros for scattering of particles with the same sign 
(i.e. e + u scattering in our case) and zeros outside the physical region for e + d scattering, 



consistent with the results derived in Refs. [ O, [12j . However we also found that for a short- 
lived resonance (Flq — > 00) and for the Standard Model there were radiation zeros also for 
e + d scattering within the physical region. Both types of Standard Model e + q — > e + q + 7 
radiation zeros will be the focus of the present study. 

From an experimental point of view the detection of photons in the final state is highly 
non-trivial. The rates are small (suppressed by 0(a) compared to the total cross section) 

1 In fact in general this is only true at tree level, see Ref. pOj. 
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and the photons must be well-separated from the beam and from the other final-state 
particles, and contained within the detector. The basic question is whether the radiation 
zeros of the scattering amplitude correspond to 'detectable' photons at HERA. In this study 
we will present results for typical values of the DIS variables y and Q 2 which correspond 
to observable quark jets and scattered positrons. For these values we will investigate the 
location of the radiation zeros for photons with an energy greater than 5 GeV. 

The paper is organised as follows. We first consider soft-photon emission and derive 
analytic solutions for the location of the radiation zeros in the eq c.m.s. frame. We then 
show how the transition from soft- to hard-photon emission shifts the position of the zeros. 
Finally we move to the HERA lab frame to see where the zeros occur in the detector. We also 
compare our exact matrix-element results with an approximate calculation in which photon 
emission is included in the colhnear approximation, which could correspond for example to 
a parton-shower implementation of such emission. This model has no radiation zeros and 
serves as a benchmark for the amplitude suppression in the exact result. 



2 Radiation zeros in eq — > eq^ scattering 

In the following we shall study the reactions 

e + (p 1 )u(p 2 ) -> e + (p 3 ) w(p 4 ) +7(fc), (1) 
e + ( Pl )d(p 2 ) -> e + (p 3 )d{p 4 ) + 1 (k). (2) 

Other scattering combinations (e + u, e~u, . . .) can be obtained from these basic processes by 
readjusting the charge factors. The expression for the matrix element squared (summed and 
averaged over spins) may for example be obtained by crossing the expression for e + e _ — > 
+ 7 given in Ref . |L6) . In terms of the four-momenta defined in Eqs. (|l]Jl) the matrix 



element for massless quarks and leptons is 

IT7 12/ + + | \ 6 2 (gl • P^) 2 + (gg • Pif + (gl • PA) 2 + (gg • Psf 

\M 3 \ (e g - e q + 7 ) = e e q j- - ^ M , (3) 



with 



\jq u = elm - e, {[12] + [34] - [14] - [23]} + [13]. (4) 
We have used here the following short-hand notation for the eikonal factors: 

= (pJkKpyky (5) 

The expression in Eq. (H) — the antenna pattern of the process — contains colhnear (k-p^ — > 
0) as well as infrared (cj 7 = E\. — > 0) singularities. It is this factor which vanishes for certain 
configurations of the momenta. Note that we only take the neutral current 7*-exchange into 
account as the antenna pattern in Eq. (f|) is independent of the exchanged particles as long 
as they do not themselves emit photons. This approximation will influence the cross section 
rate slightly at high Q 2 , but will not affect the position of the radiation zeros. 
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2.1 Type 1 radiation zeros 

To see under what conditions JFg M vanishes, we first recall the 'single-photon theorem' from 
Ref. JT3 which states that the amplitude vanishes when the charge-weighted scalar products 
Qi/{Pi ' k) are equal. If we denote the common value by A, then 

[ij] = {QiQ j )- 1 >?Pi-p j (6) 

and it is straightforward to show by substitution in Eq. (H) that this gives jFg M = 0. In the 
present context, the equality of the charge- weighted scalar products corresponds to 



Cq 1 Gq 



Pi ■ k p 2 ■ k P3 ■ k Pi - k 



(7) 



We can obtain a simple analytic solution to these equations by taking the soft-photon limit 
in which uo^/Ei — > 0. In this limit we have simple two-body kinematics for the quarks and 
leptons, pi +p 2 = P3 +Pa- If we work in the e + q c.m.s. frame, and define 9 2 , #4 to be the angle 
between the photon and the incoming and outgoing quarks respectively, then the equations 
(0) become 

1 e„ 1 e„ 



1 + z 2 1 - z 2 1 + z A 1 - Z4 
where Zj = cos 6^. Equivalently, 

z 2 = z 4 = i^i (9) 
1 + e q 

A necessary condition for such a solution to physically exist is e q > (=>■ \zi\ < 1), i.e. e + u 
or e + d scattering. This reproduces the well-known result for scattering of particles with 



the same sign of electric charge, as discussed in Refs. p0| . We call these Type 1 radiation 
zeros. By itself, however, the condition e q > is not sufficient to guarantee a zero in the 
scattering amplitude. The equation z 2 = z^ can only be satisfied for certain configurations 
of the final-state particles. To see this, we introduce an explicit representation of the c.m.s. 
four-momenta: 

(10) 

(11) 
(12) 

(13) 

cosfl 7 ). (14) 

These variables are illustrated in Fig. ffl. It is straightforward to show that the conditions 
for JF^ M = defined in Eq. (£5|) correspond to 

cos l-\^l (15) 
7 1 + e' 



Pi = 


2 


(1,0,0,-1) 




P2 = 


2 


(1,0,0,1), 




pl = 


v§ 


(l,sin6 9 ,0, 


cos Q q ) , 


P3 = 


2 


(1, -sin6g, 


0, — COS Q q ] 


P = 


^7 ( 


l,sin 9 1 cos^ 


4 7 , sin 6^ sin 
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and 



± arccos 



Aan(8 q /2) ' 
\ tan # 7 



78.46° and for eg = +1/3 at 9 1 



(16) 



60 c 



Thus for e u = +2/3 we find radiation zeros at Q 1 
We present the positions of the radiation zeros (</> 7 ,0 7 ) for process ([!]) (e + u scattering) 
in Fig. Qa. Note that the requirement of a physical solution for 7 places restrictions on 
Q q . There are two radiation zeros in the (0 7 , 6* 7 ) plane for Q q < 26 1 ~ 156.94°. The 
cones around the incoming and outgoing quarks defined by £2,-24 = 1/5 have two lines of 
intersection along which there is completely destructive interference of the radiation. Note 
also that at Q q = 2# 7 = Q q nt the radiation zeros degenerate to a single line (i.e. single point 
in (0 7 , # 7 ) space) located in the scattering plane (0 7 = 0°). There are no radiation zeros for 
Q q > 26 1 ~ 156.94°. Finally, for Q q = 0° there is an infinite number of radiation zeros ('null 
zone') located on a cone around the beam line with opening angle # 7 . 



2.2 Type 2 radiation zeros 

The processes (0,0) exhibit a second class of radiation zeros, which we call Type 2, which 
do not satisfy the 'single-photon theorem'. These zeros are located in the scattering plane at 
7 = 0° and 7 = 180°. The corresponding 9 1 values may be calculated straightforwardly in 
the soft-photon approximation as a function of the quark charge e q and the quark scattering 
angle Q q . The result is 



cos 6~ 



1 (l - efj (1 + cos0y) ± y / A 7 (e 5 ,cos Q q 

2 (l~e q ) 2 



(17) 



with 



A 7 (e ff , cos Q q ) = [(e 2 q - l) (1 + cos 6,)] 2 - 4 (1 - e q f (e 2 q cos Q q + 2e q + cos Q q ) . (18) 

The condition A 7 (e q ,cosQ q ) > constrains the range of e q for which physical zeros exist. 
In terms of the polar angle Q q we have 



cos Q q + 3 - 2^2(1 + cos Q q ) 
00 < e q < y -— < 1, 

1 — COS Hg 



or 



1 < 



cos 6 o + 3 + 2J2(l + cos0 



cos 8, 



< e„ < +00, 



(19) 



(20) 



the latter being actually redundant since Standard Model quarks have \e q \ < +2/3. From 
Eq. ([19"D we obtain constraints on the quark scattering angle Q q for particular flavours of 
quark. There are radiation zeros for all e q < and for positively charged quarks in a limited 
range of Q q . We summarise the results in Table |1|. Note that e + u scattering has both Type 
1 and 2 zeros. However, the latter are located very close to the beam direction, making their 
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e + u 


e q = +2/3 


cos & q < 7r — arccos (||) 


e> q > 157° 


e + d 


e q = +1/3 


cos Q q < —\ 


e q > i2o° 


e + d 


e q = -1/3 


V COSOg 


ve, 


e + u 


e q = -2/3 


V COSOg 





Table 1: Ranges of the quark scattering angle Q q , for different quark charges, for which 
radiation zeros exist. Note that for e q < there are always two radiation zeros in the 
scattering plane for 7 = (0°, 180°) with the 9 1 value given by Eq. (17). 



observation difficult in practice. They also require very high Q 2 (back-scattered quarks) 
and therefore have a small event rate. The positions of the Type 2 zeros for e + d scattering 
are shown in Fig. ^|b as a function of Q q . Finally, Table |2| lists the numerical values of the 
radiation zero angles (0 7 ,6* 7 ) for several values of Q q . 



0, 


e + d scattering 


e + u scattering 


30° 


(0°, 76.12°) 


(180°, 46.12°) 


(-86.86°, 78.46°) 


(86.86°, 78.46°) 


45° 


(0°, 84.98°) 


(180°, 39.98°) 


(-83.23°, 78.46°) 


(83.23°, 78.46°) 


90° 


(0°, 114.29°) 


(180°, 24.29°) 


(-78.22°, 78.46°) 


(78.22°, 78.46°) 



Table 2: Position of the radiation zeros (0 7 , # 7 ) for three different quark scattering angles 
Q q , in the soft-photon approximation. 



2.3 Radiation zeros for arbitrary photon energies 

The analytic results obtained above use the soft-photon approximation. However radiation 
zeros of both types exist for all photon energies and can be located using numerical tech- 
niques. We continue to work in the e + q c.m.s. frame but now use exact 2 — > 3 kinematics. 
Without any essential loss of generality, we can keep the direction (Q q ) and the energy (E' q ) 
of the outgoing quark fixed and vary the direction and energy of the outgoing photon, con- 
structing simultaneously the four-momentum of the outgoing positron to conserve energy 
and momentum. The new four-vectors of the outgoing quark, lepton and photon momenta 
are then 

v l = (i, an e ff ,o, cose,), (21) 

P"s = pf+rf-rf-A^, (22) 
= oo 7 (1, sin^ 7 cos0 7 ,sin0 7 sin0 7 ,cos^ 7 ) . (23) 



5 



Once again we obtain a vanishing matrix element in Eq. ([|) if the antenna pattern JF^ M 
of Eq. (f|) is zero. For Type 1 radiation zeros, the single-photon theorem again leads to the 
conditions in Eq. (|?p. The equality of p% ■ k and P2 ■ k leads immediately to Eq. ([T5|), i.e. 
the radiation zeros are at fixed # 7 independent of the photon energy. However the azimuthal 
angle 7 does vary with cj 7 , since the supplementary condition z<i = 24 only applies in the 
cij 7 — > limit. For Type 2 zeros, it can be shown that the condition 7 = 0°, 180° again 
applies for arbitrary ou 1 , i.e. the zeros are always located in the scattering plane. 

In Fig. |3| we show the dimensionless quantity A/"g M = w 7 JFg M for different photon energies 
and fixed final-state quark kinematics. The figures (a) and (b) correspond respectively to 
slices through the (0 7 , # 7 ) plane according to the positions of the soft-photon Type 1 and 2 
radiation zeros of the previous sections. As the photon energy increases, there is a systematic 
shift in the positions of the zeros. As radiation zeros are semi-classical effects due to de- 
structive interference, it is easy to understand that fixing the position of the outgoing quark 
and simultaneously increasing uj 1 shifts the interference regions between the participating 
charged particles as the outgoing positron must balance energy and momentum and thus 
changes its relative orientation. Thus the asymmetric uj 1 dependence of the two radiation 
zeros in Fig. ^a is due simply to our choice of fixing the final-state quark direction rather 
than the direction of the scattered positron. The zero in the quadrant between the (fixed) 
incoming positron and outgoing quark directions is relatively insensitive to the changes in 
the positron direction induced by varying cu 7 . The other zero follows the direction of the 
outgoing positron as u; 7 increases. The same effect also explains the symmetric dependence 
of the two radiation zeros for the process e + u — > e + u + 7. The zeros are located symmet- 
rically above and below the scattering plane and are influenced equally by changes in the 
scattered positron direction. 

Figs. ^a,b show the positions^ of the radiation zeros Z d and Z u for the two processes as 
a function of the photon energy at various fixed 9,j. Z d (Q d , ^7) is located in the quadrant 
between the outgoing positron and the incoming d quark and Z u (Q u ,u^) is located in the 
quadrant between the outgoing u quark and the outgoing positron (see Figs. ^|a,b). The 
values on the axes at uo^ = coincide with the analytic results obtained previously (see 
Table |). 

The dashed lines in Figs. |]a,b are simple polynomial fits. For e + d scattering we fit 9 1 for 
fixed 7 = 180° using a quadratic polynomial, 

Z d (Q d , cu 7 ) = Z° d (Q d ) + d l u 1 + d 2 u 2 v (24) 

where Z d (Q d ) corresponds to the soft-photon results listed in Table |2|. The radiation zeros 
for e + u scattering (i.e. a fit for 7 at fixed 9 1 = 78.46°) can be approximated by a first-order 
polynomial 

Z U (Q U , w 7 ) = Z®(Q U ) + «i^ 7 . (25) 

The results of the fit are presented in Table ^. 

As a final exercise in our c.m.s. studies we calculate the differential cross section for the 
two subprocesses. The general form of the differential subprocess cross section in the e + q 

2 The exact locations of the zeros are determined by a numerical procedure. 
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d-quarks 


u-quarks 




di 


d 2 


x 2 




x 2 


30° 


0.576 


0.015 


5.91 


0.59 


5.3 x 10" 2 


45° 


0.012 


0.014 


0.76 


0.40 


3.0 x 10~ 3 


60° 


0.149 


0.002 


0.01 


0.28 


1.3 x 10~ 4 



Table 3: Fits for the uj 1 dependence of the two selected radiation zeros shown in Figs. $a,b 
according to the definitions given in Eqs. ( f^ , \2dj ). 



c.m.s. frame may be written as 

— — — eg -> eg + 7 = — — / du y 9 7 M 3 2 (eq -> eg + 7), 26 



jCUt 



where 



e; = ^ g ~ 2 ^ , ~ . (27) 



x s — 2a; 7 (l — cos 97; 



The integration over lo 1 smears out the radiation zeros to form a sharp dip in the cross 
section. Since the cross section decreases rapidly with increasing cj 7 , the dip is close to 
the location of the zero corresponding to fixed u; 7 = ^ ut . The distributions for the two 
subprocesses are shown in Figs. |5|a,b for oo^ ut = 5 GeV. Note that we have also imposed an 
angular cut around the beam line of 5° in Fig. [|a. The transition from radiation zeros to 
radiation dips can be seen by comparing Figs. | and || Choosing larger values of ti; 7 ut shifts 
the radiation dips to higher values of 7 and 9 1 at the same time decreasing the overall value 
of the subprocess cross section. 



3 Radiation zeros at HERA 

In this section we shall discuss the possible observation of radiation zeros at HERA. To do 
this we modify the previous calculation by (a) moving to the HERA lab frame, (b) including 
the parton distribution functions, and (c) summing over all flavours of quarks in the initial 
state. 

In neutral current DIS the cleanest way to reconstruct the kinematics of a given event 
is by measuring the energy E' e and the laboratory angle Gg ab of the outgoing positron. In 
terms of the Bjorken scaling variables x and y we may write (see for example Refs. [17], |T 

E: 



2E, 



(l-cosGf 3 ), (28) 
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* = ^|(l + cosef), (29) 

Q 2 = xys, (30) 

where E p is the energy of the incoming proton and s = 4E e E p is the c.m.s. energy of the e + p 
system. The polar angle of the positron 0|f b is defined with respect to the incident proton 
beam direction. The precision of the y measurement typically degrades as 1/y, and thus one 
naturally assumes y > 0.05 (18[ . 



Since we are interested in DIS events with an additional hard photon emitted at different 
angles in phase space, the natural quantity to consider is the triple-differential cross section 
d 3 a/dydQ 2 d^ h . In the HERA lab frame this is given by 

d 3 cr . 1 

(e + p - e + q + 1 + X) = — — -£ 



dydQWf * * 1 ' 256vr 4 s 

/UJ 
dUl £ (oy x -2v-k) ^ 3|2(e9 eq + 7) q2); (31) 



^cut 



where 



£ q (Q 2 /x-2p.k) 
Q 2 - 2q ■ k 



^ q Q 2 /x-2p-k~ X - (32) 
In the calculations which follow we choose E e = 27.5 GeV, E p = 820 GeV and neglect all 



y 


X 


Cjlab 

e 


K 


0.20 


0.55 


52.4° 


112.9 GeV 


0.40 


0.27 


46.2° 


107.4 GeV 


0.60 


0.18 


38.4° 


101.9 GeV 


0.80 


0.14 


27.6° 


96.4 GeV 



Table 4: Typical values of the scattered positron energy and angle for our parameter choice 
Q 2 = 10 4 GeV 2 and different values of y. 



quark and lepton masses. We again take uo™ = 5 GeV for the lower limit of the photon 



energy. For the quark distribution functions f q / p (£ q ,Q ) we use the MRS(A') set of partons 



introduced in Ref. [0, with QCD scale parameter A^- -4 = 231 MeV corresponding to 
a s (M|) = 0.113. In order to stay in the valence-quark scattering region (i.e. large £ q ), 
where we expect the radiation zeros to be most visible, we choose Q 2 = 10 4 GeV 2 and 
y G [0.1,1.0]. Typical values for x and the positron variables Gg ab and E' e are listed in 
Table |. 

As we move from the e + q c.m.s. frame to the HERA lab frame, all four-momenta are 
boosted along the beam direction. Although this has no effect on the azimuthal angles, the 
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polar angles and hence the locations of all radiation zeros, in particular # 7 , are changed. The 
simplest consequence of this is that the e + d radiation scattering zeros remain located in the 
scattering plane at 7 = 0° and 180°. To find the locations of the radiation zeros for process 
(§) we therefore fix 7 = 0° and numerically determine their positions in 6 I 7 . 

3.1 Radiation zeros for d quark scattering 

In Fig. |^ we present the cross section of Eq. ([H]) for the process e + + p — > e + + jet + 7 + X 
via e + d — > e + d + 7 scattering (dashed line) as well as via the sum over all subprocesses 
e + q — > e + q + 7 with q = u,d, s, u, d and s. We have chosen to focus on the radiation zero 
located between the incoming quark and outgoing positron. We fix Q 2 = 10 4 GeV 2 and vary 
y from y = 0.2 in Fig. §a to y — 0.8 in Fig. |6|d, which corresponds to x values in the region 
0.1 < x < 0.6 (cf. Table |J). Again we observe radiation dips instead of radiation zeros due 
to the integration over the photon energy. Increasing y pulls the radiation dips closer to 
the beam line and thus makes their observation more difficult. Already at y = 0.2 the e + d 
radiation dip in Fig. ^|a is only about 14° (cos# 7 ~ 0.97) from the beam line, and gets even 
closer with increasing y. Note that we impose a cut of 5° around the beam line. Increasing y 
means decreasing the polar angle of the outgoing positron Bg ab (cf. Table |J). Thus the zone 
of destructive interference approaches the beam line as the e + approaches the beam line. 
The conclusion is that observation of the radiation dips in the sector between the incoming 
quark and outgoing e + in high-Q 2 events is only possible for small values of y. 

The second radiation zero we found in our studies was located between the incoming 
positron and the outgoing quark. In Fig. [7| we display this region again for processes only 
involving d quarks (dashed lines) as well as for processes involving all light quark and an- 
tiquark flavours. The obvious singularities in Figs. [7|a-d are caused by collinearity of the 
photon with the outgoing quark. Now the problem is that the zeros are close (always within 
10°) to the outgoing quark jet, even though the radiation dips here are well separated from 
the beam line (~ 35° for y = 0.6). 

A more serious problem evident in Figs. |6] and [7] is the enormous background from the 
other quark scattering subprocesses, which completely fills in the radiation dip. We observe 
a ratio (away from the singularities) of signal/background ~ 1/ (200 — 300). The dominance 
of the u quark contribution is striking. For the given values of x and thus £ g (cf. Eq. (p2|)) 
we find the following ratios for the MRS (A') parton distributions at Q 2 = 10 4 GeV 2 : 

£ q = 0.1 : - u(Q : d(Q : d(Q : u(Q ~ 100 : 60 : 22 : 15, (33) 

^ = 0.6 : - u{Q : d(Q : dfe) : u{Q ~ 100 : 17 : 1 : 1. (34) 

In addition to these parton distribution factors there are the usual quark charge squared 
(e 2 ) factors from the leading order eq —>■ eq scattering, which further enhance the u-quark 
contribution. Note that the s-quark contribution plays a minor role; it is roughly 70% of 
the u contribution at £ g = 0.1 and comparable to the latter at higher values of £ 9 . Even 
though d, s and u quarks all yield radiation dips in the scattering plane (the d- and s-quark 
zeros coincide) none of these are likely to be observable. The only possibility might be to try 
to flavour-tag the d or s quark jets, for example by selecting only those jets with a leading 
negatively charged track. 
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3.2 Radiation zeros for u quark scattering 

According to the parton distribution hierarchy presented in the previous section we might 
expect that the Type 1 radiation zeros, which we identified with the traditional radiation 
zeros already discussed in the literature, are the most promising for detection. We recall 
that in the soft-photon limit and in the c.m.s. frame these zeros are located at fixed polar 
angle cos# 7 = 1/5 (cf. Eq. (|1|)). Their position in 7 may then be directly computed 
using Eq. (|i"6|) . We found that they are located well outside the scattering plane (except for 
Q q = 29 y = 2cos~ 1 (l/5)) as discussed earlier. Integrating over the photon energy uj 1 and 
using exact 2 3 kinematics slightly shifts the position of the corresponding radiation dips. 
The uj 1 dependence for different kinematical situations was shown in Fig. |4]b. 

Moving to the HERA lab frame boosts the polar angles and changes the position of 
the radiation dips for e + u — > e + u + 7. In Figs. ||a c we show the differential cross section 
of Eq. (|3~ID for this process over the full (0 7 , # 7 ) space. As before we fix Q 2 = 10 4 GeV 2 
and chose the three y values: 0.2, 0.4 and 0.6. We impose cuts of 5° around the beam 
line (by definition located at Q 1 = 0° and 180°) and cut the differential cross section at 
da < 10~ 4 pb/GeV 2 to avoid the collinear singularities along the directions of the outgoing 
e + (located at 7 = 0°) and the outgoing u quark (at <p 1 = ±180°). We see that the positions 
of the zeros are still symmetric in 7 , as expected. Note that since the collinear singularities 
and the radiation dips tend to concentrate around small values of # 7 , we have introduced a 
logarithmic scale for # 7 in the three-dimensional plots of Fig. || 

We can numerically locate the positions of the radiation dips in the (0 7 , # 7 ) phase space 
for our different choices of y: 

y = 0.2 -> 7 ^ ±97.2°, 6 1 ~ 20.6°; 
y = 0A -> 7 ~ ±100.4°, 9 7 ~ 24.9°; 
y = 0.6 -> 7 ~ ±102.5°, 6 1 ~ 25.2°. (35) 

It is straightforward to verify that the radiation dips, if projected onto the scattering 
plane, lie within the quadrants between the incoming (outgoing) e + and the outgoing (in- 
coming) quark, the zone of destructive interference. As Fig. ||] shows, the radiation dips are 
clustered quite close to the (beam) direction of the incoming quark (# 7 = 0°) which is par- 
ticularly true for high-Q 2 events (back-scattered positron). As we have already pointed out, 
they are also within 10° (in # 7 ) of the final-state quark jet. However, they are well-separated 
from the outgoing particles when the 7 angle is taken into account. It will be very important 
to perform realistic simulations of these photon radiation events, including jet fragmentation 
and detector effects, to see whether the dips are indeed observable in practice. 

Finally, in Fig. || we show the 7 dependence for slices through the # 7 values of Eq. (|35| ) 
which define the numerical location of the radiation dips of Fig. We show the contributions 
of u quarks only, as well as the contributions from all light flavours (i.e. u, d and s quarks 
and antiquarks). At the critical values of 7 (again given in Eq. fl3"5|)) the obvious dips for 
pure u-quark scattering are somewhat filled in by the other 'background' (mainly (i-quark) 
processes — the cross section at the bottom of the dip is increased by about two orders of 
magnitude — although they are still significant. 
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3.3 Radiation zeros and 'parton shower' models 



To gauge the quantitative significance of the radiation zeros described in the previous sec- 
tions, and in particular to factor out the effects of phase space constraints on the distribu- 
tions, it is useful to make comparison with an approximate calculation in which radiation 
zeros are absent. Parton shower Monte Carlo programs, such as Herwig [^(J or Pythia 
PH , are based on the principle of the leading-pole (collinear) approximation. In particular 



they do not usually include the interference effects which are crucial for producing radiation 
zeros in the scattering amplitudes. We can easily emulate such models by removing the 
interference terms from the antenna pattern in Eq. (i.e. the terms linear in e q ): 

= ej[24] + [13], (36) 

The approximate matrix element thus obtained still contains the correct leading collinear 
singularities when the photon is emitted parallel to the incoming and outgoing quarks and 



leptons. In Fig. 10 we present the ratio 



approx 

RU , = ( e+u "> e+u + , (37) 



where d 3 a a,ppTOX /dydQ 2 dQ l ^ b includes the antenna pattern without interference terms, as 
defined in Eq. (|3~6p. Again we slice through 7 at the values 9 1 of Eq. (|3~5| ) where we 
numerically located the positions of the radiation dips for each y value. Note that away 



from the dips the ratio is 0(1), as expected. However Fig. 10 also shows that close to 
the dips the approximate cross section is up to three orders of magnitude larger than the 
exact result, for all y values. In these particular regions of phase space, therefore, such 
'parton-shower' models would dramatically overestimate the photon emission cross section. 



4 Conclusions 

The scattering amplitude for the process eq — > qe + 7 vanishes for certain configurations 
of the final-state momenta. In this paper we have studied these radiation zeros and in 
particular their observability at HERA. In addition to the well-known class of (Type 1) 
same-charge zeros, which have been discussed in the pioneering work of Refs. [|ll|] and [|Tj| , 
we have discovered a second class of (Type 2) zeros located in the eq scattering plane. We 
have so far been unable to find a theorem which leads to conditions for the existence of such 
zeros in more general scattering processes. 

Experimentally, one might hope to be able to measure the four-momenta of the final- 
state lepton, quark (jet) and photon sufficiently accurately that the kinematic configurations 
which lead to zeros could be reconstructed. However a more realistic approach, which we 
have adopted here, is to study DIS + photon events for fixed lepton variables y and Q 2 
and for a range of photon energies above a given threshold. This leads to sharp radiation 
dips instead of zeros. We performed such a study using the HERA lab frame. Although 
the radiation dips, i.e. the photon directions for which the cross section has a minimum, of 
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both types are quite well separated from the beam direction and from the final-state jet, 
the e + d scattering dips are completely swamped by the contributions from the other quark 
scattering processes. The e + u (Type 1) dips offer a more promising hope of detection, since 
e + u scattering is the dominant subprocess at high x. With sufficient statistics, we would 
expect such dips in the cross section to be observable at HERA. 
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Figure 1: Parametrisation of the kinematics for e 4 '{pi)q{p2) — > e + (j»3)g(j»4) +7(fc) scattering 
in the e + q c.m.s. frame. The orientation of the photon relative to the scattering plane is 
denoted by # 7 and <p T 
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Figure 2: The position of the radiation zeros as a function of the quark scattering angle <d q 
for soft-photon emission in (a) e + u — > e + u + 7 and (b) e + d — > e + d + 7 in the (0 7 , 9 y ) c.m.s. 
phase space of the soft photon. 
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Figure 3: The dimensionless antenna pattern Af$ M = u}^J-g M for (a) e + d — > e + d + 7 (at 

fixed 7 = 0°) and (b) e + u — > e + u + 7 (at fixed # 7 ~ 78.46°) for different photon energies 
(ou 1 = 10,20,30,40 GeV). The outgoing quark direction is fixed at Q q = 90° with energy 
E' q = 100 GeV. The directions of the incoming and outgoing quarks and leptons are indicated. 
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Figure 4: The positions of the radiation zeros Z q as a function of the quark scattering angle 
G g and the photon energy u; 7 for (a) e + d scattering and (b) e + u scattering. The analytic 
results for the soft-photon limit {u 1 — > 0) are summarised in Table [|. The dashed lines are a 
polynomial fit in the photon energy for given Q q . Note that in (a) we employ a second-order 
fit whereas in (b) a first-order fit is sufficient. The fit parameters are listed in the text. 
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Figure 5: The subprocess differential cross section d 2 a/dQ^dQ q of Eq. (|26| ) for c.m.s. e + d 
scattering (solid lines) and e + u scattering (dashed lines). We again choose those slices 
through the photon parameter space (0 7 , # 7 ) that contain radiation zeros in the soft limit 
(i.e. a choice of 7 = 0° in (a) and cos# 7 = 1/5 in (b)). Note that we integrate over 
the photon energy cj 7 and fix the position of the outgoing quark at Q q = 90° with energy 
E' = 100 GeV. In (a) we impose an additional angular cut of 5° around the beam line. 



19 





Figure 6: The positions of the radiation dips for the process e + d —>■ e + d + 7 (dashed lines) 
for different values of y for 6* 7 G [5°, 60°] (hemisphere of outgoing e + ). A cut of 5° around 
the beam line is imposed. The radiation zeros and thus the radiation dips for this process 
are again located within the plane (0 7 = 0°). The Q 2 value is 10 4 GeV 2 . The solid lines 
show the sum of the contributions from u, d, s quarks and antiquarks. The divergences in 
the plots show the positions of the outgoing e + at cos 9 1 = cos 6g ab with the values for 6jf b 
given in Table |j. 
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Figure 7: The positions of the radiation dips for the process e + d — > e + d + 7 (dashed lines) 
for different values of y for # 7 e [5°, 60°] (the hemisphere of the outgoing quark). A cut of 
5° around the beam line is imposed. Note that 7 = 180°. 



21 




Figure 8: The differential cross section of Eq. (|31| ) for the process e + u — > e + u + 7 in the 
(0 7 , 9ry) phase space of the emitted photon. The corresponding contour plots are shown on 
the right-hand side. For fixed Q 2 = 10 4 GeV 2 we vary y (defined in Eq. (|2"ED) from y = 0.2 
in (a) and y = 0.4 in (b) to y = 0.6 in (c). The kinematic variables x, E' e and 0g ab for each 
y value can be read off from Table |j. Note that we introduce in the surface plots on the 
left-hand side a logarithmic scale in 6* 7 . The radiation dips are symmetric in 7 . Again we 
impose a 5° cut around the beam line, and thus 6* 7 G [5°, 175°]. 
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Figure 9: The differential cross section of Eq. fl3~T|) for three different y values at the position 
of the radiation dips # 7 shown in Fig. |8] as a function of the azimuthal angle 7 . We show 
the process e + u — > e + u + 7 (dashed lines) as well as the contribution (solid lines) of all light 
quark flavours (u,d,s quarks and antiquarks). 
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Figure 10: Same as Fig. || but now for the process e + u — > e + u + 7 only. is the ratio of 



the differential cross sections of Eq. (|31|) without and with interference terms, see Eq. 
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